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Abstract. In this paper we analyze so-called Parisian ruin probability 
that happens when surplus process stays below zero longer than fixed 
amount of time £ > 0. We focus on general spectrally negative Levy 
insurance risk process. For this class of processes we identify expression 
for ruin probability in terms of some other quantities that could be pos- 
sibly calculated explicitly in many models. We find its Cramer-type and 
convolution-equivalent asymptotics when reserves tends to infinity. Fi- 
nally, we analyze few explicit examples. 
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1 Introduction 

In risk theory we usually consider classical Cramer-Lundberg risk process: 



is a compound Poisson process. We assume that Ui, (i = 1,2, ...) are i.i.d dis- 
tributed claims (with the distribution function F and tail F := 1 — F ). Arrival 
process is a homogeneous Poisson process Nt with intensity A. Premium income 
is modeled by a constant premium density p and the net profit condition is then 
\v jp < 1, where E(C/i) = v < oo. Lately there has been considered more gen- 
eral setting of a spectrally negative Levy process. That is, X — {X t }t>o is a 
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X t — x + pt — S t 



(1) 



where x > denotes the initial reserve, 




(2) 



process with stationary and independent increments with only negative jumps. 
We will assume that process starts from X n = x and later we will use conven- 
tion P(-|Xo = x) = P x ( - ) and Po = P. Such process takes into account not only 
large claims compensated by a steady income at rate p > 0, but also small per- 
turbations coming from gaussian component and possibly additionally (when 
IIx(— oo,0) = 00 for jump measure Tlx of X) compensated countable infinite 
number of small claims arriving over each finite time horizon (see e.g. |15)V 

One of the most important characteristics in risk theory is a ruin probability 
defined by P(tJ" < 00) for = inf{i > : X t < 0}. In this paper we 
extend this notion to called Parisian ruin probability, that occurs if process X 
stays below zero for longer period than a fixed £ > 0. Formally, we define the 
excursion below zero: 

ft = sup{s < t : l(x s >0)!(x t <0) = !}• 

Parisian time of ruin is given by 

r c = inf{i > : t - ft > C} 

and Parisian ruin probability we define as: 

P(t c < oo\X =x)= P x (t c < 00). 

The case Q = corresponds to classical ruin problem and we do not deal with 
this case here. The name for this problem comes from Parisian option that 
prices are activated or canceled depending on type of option if underlying asset 
stays above or below barrier long enough in a row. We believe that Parisian ruin 
probability could be better measure of risk in many situations giving possibility 
for insurance company to get solvency. Parisian ruin probability was already 
considered by Dassios and Wu [7] where Parisian ruin probability was found 
for classical risk process <m with exponential claims and for Brownian motion 
with drift. In Dassios and Wu [8] Cramer- type asymptotics was found for (Q}. 
In this paper using the fluctuation theory we show that these results could be 
extended to the case of a general spectrally negative Levy process. In particular 
we show how to deal with risk process being sum of independent classical risk 
process (JTJ) and Brownian or a-stable motion perturbations. Additionally, we 
derive asymptotics of Parisian ruin probability in convolution-equivalent case, 
that is when claim size has a heavy-tailed distribution. 

This paper is organized as follows. In Section [5] we introduce basic notions 
and notations. In Section [3] we give main representation of Parisian ruin prob- 
ability. In Sections 2] and [5] we give asymptotics of Parisian ruin probability in 
Cramer and convolutions-equivalent case, respectively. Finally, in Section [6] we 
analyze few particular examples. 

2 Preliminaries 

In this paper we consider a spectrally negative Levy process X = {X t } t >o, that 
is a Levy process with the Levy measure ILjf satisfying nx(0, 00) = 0. Assume 
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that Xq — x > and that X t — > oo as t — > oo a.s. (that is, reserves of the 
insurance company increases to infinity a.s.). This assumption excludes also 
case of a compound Poisson process with negative jumps. We assume also that 
KXi < oo. Note that EXi > by drift assumption. With X we associate 
the Laplace exponent <p(/3) ■= j\ogK(e l3Xt ) defined for all /3 > and function 
= sup{6> > : tp(9) = q}. We will consider also dual process X t = —X t 
which is a spectrally positive Levy process with jump measure II o (0,y) = 
IIx (— y,0). Characteristics of X will be indicated by using a hat over the 
existing notation for characteristics of X. For process X we define ascending 
ladder height process (L _1 ,iJ) = {(L^~ , H t )}t>a- 



Lt 



inf{s > : L s > t} if t < 

oo otherwise 



and 

H t := 



X j -i if t < Loo 
oo otherwise, 



where L — {L t } t >o is a local time at the maximum (see [T71 p. 140]). Recall 
that (L~ l ,H) is a bivariate subordinator with the Laplace exponent k(6,(3) = 

— j logE (e~ 9L t ~^ Ht l{t<L 00 }J and with the jump measure IT#. We define 

descending ladder height process (L _1 ,1Z") = {(L^ 1 , H t )}t>o with the Laplace 
exponent /?(#,/?) constructed from dual process X. Recall that has expo- 
nential distribution with parameter k(0,0). Moreover, from the Wiencr-Hopf 
factorization we have: 

see [HI p. 169-170]. Hence 

£(0,0) = ^(0+). (4) 
We introduce a potential measure U defined by 



U{dx,ds)=\ W(Lt l eds,H t £ 
Jo 



dx)dt 



with the Laplace transform JL > 2 e 6s ^ x U(dx,ds) = 1/k(8,/3) and renewal 
function: 

U(dx) = f U(dx, ds) = E ( [ l {Htedx} dt) . 

J[0,oo) \J0 / 

For spectrally negative Levy process upward ladder height process is a linear 
drift and hence the renewal measure is just Lebesgue measure: 

U{dx) = dx. (5) 
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Moreover, from (J3|), 



"U(dz) = Wy (6) 



see [13 p. 195]. 

We will also use the first passage times: 

t~ = mf{< > : X t < x}, r+ = inf{< > : X t > x}. 

Suppose now probabilities {P^jxeR corresponds to the conditional version of P 
where Xq = x is given. We simply write Po = P. We define also Girsanov-type 
change of measure via: 

= £M (7 ) 

for any c for which Ke cXl < oo, where St (c) = exp{cX t — ip (c) t} is exponential 
martingale under P x and Tt is a natural filtration of X. It is easy to check that 
under this change of measure, X remains within the class of spectrally negative 
processes and the Laplace exponent of X under P° is given by 

<p c (9) = <p(6 + c)-cp (c) (8) 

for 9 > — c. Similarly, by © on P c exponent of descending ladder process equals, 

«? fn R\ ^ (Z 3 + c ) - y (c) /q\ 

KC{0,P) = ^ (9) 

since on an account of positive drift of X we have $(0) = 0. 

For q > 0, there exists a function : [0, oo) — > [0, oo), called q-scale 

function, that is continuous and increasing with the Laplace transform 



e -6x 



WW(y)dy = (<p(d)-q)-\ 6 > $(q). (10) 



We denote W^°'{x) — W(x). Domain of is extended to the entire real axis 
by setting W^(y) — for y < 0. For each y > 0, function q — > W^(y) may 
be analytically extended to q e C. Moreover, let 

rv 

ZW(y) = l + q W [q \z)dz. 
Jo 

It is known that: 

E x (e- 6T v,T+<oc) = e -*W^-), (11) 
M 6 ^V<~)) = ^(x)-^W^(x), (12) 
= + ^g^ W^^s)) (13) 



e 
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where W^ v{v)) and Z { ~ v{v)) are scale functions calculated with respect to the 
measure W . We understand 6*/$ (6) in the limiting sense for 9 — 0, so that 

v x (t q <ooj-| l if^'(o)<0. li4j 

For details see [17] or [16l Remark 3]. If F X (X T - = 0, t ~ < oo) > then we say 
about downward creeping. This is possible only when downward ladder process 
is of unbounded variation with strictly positive drift d. In this case the renewal 
function U has a strictly positive and continuous density u on (0, oo) satisfying: 

V X (X T - = 0, to < oo) = du(aO, (15) 
where u(0+) = 4; see [T7| Theorem 5.9 and Problem 5.5]. 

3 Main representation 

The main representation is given in the next theorem. 

Theorem 1 Parisian ruin probability for a spectrally negative Levy risk process 
equals: 

¥ x (t c < oo) = ¥ x (tq < oo)P(t c < oo) (16) 

+ (1 - P(r« < oo)) / P(r+ > OMto < oo, -X r - E dz) 
Jo 



P x {tq < oo) = «(0, 0)U(x, oo) = 1 - tp'(0+)W(x), (17) 

p OO />00 

/ e- 0s ds / P(t+ > s)F x (tq < oo,-X T - € dz) (18) 
Jo Jo 

= g(o 1 o)g(g 1 oo) _ «(o,o) c$Wa! r 00 e _ $(e)2/ (19) 



6» $(6>) 

Proof The equation (1161) follows from the strong Markov property and from fact 
that X has only negative jumps. Moreover, identity 

¥ x (t - < oo) = P(?+ < oo) (21) 

and observation that for x > 0: 

P(r+ < oo) = /t(0,0)*7(>,oo) (22) 
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complete proof of (JT7J) in view of (0 (for see [TO p. 187]). Finally, note 
that from (fTT|) : 

f 6 " SP(T * + >S ) dS =K 1_E ( 6 " +1 K<oo))) = J (l - • 

(23) 

Further, in case of downward creeping, distribution P^(r < oo,— X - € dz) 
might have an atom at but it does not have influence on the rhs of (|2T))) : 

e~ 0s ds I P(r+ > s)F x (tq < oo, -X T - G dz) 
Jo 

= - Q £ (l- e"*W*) P x (r - < oo, -* T - G dz) 

(24) 



; (r - < oo) _ 1 



<S>(6)X _ 

e T ° , t < oo 



Identity (|T3| completes proof of ([20)l . The equation (fl"9)) follows from following 
equality: 

V X (-X T - G dz, r ~ < oo) = P(^ ? + - x G dz, ?+ < oo) (25) 
and [Til Problem 5.5]. 

□ 

Remark 1 Note that 



( X ) dx = V {p +mr\ (*) = 1 - y m d y . 

Remark 2 Consider now particular case of spectrally negative Levy process: 

X t =x+pt-S t +aB t , (26) 

where p > 0, a > 0, St is a subordinator of bounded variation with jump 
measure 11^ and Bt is a Brownian motion independent of St- Moreover, since 
St has bounded variation, then: 

/O />oc 
zHx(dz) = / U -g(z) dz < oo, 
-oo JO 

where IL^ (z) = II^(z,oo). Note that p = ES\/p < 1 because EXi > 0. We 
recall now Pollaczek-Khintchine formula (see (T3l Theorem 3.1]): 

oo 

F x (to < oo) = k(0,0)U(x,oo) = 1-(1- Po )J2po (V (n+1) * * M n *) (x), (27) 

n=0 

where 

M{dz) = —U s (z)dz, (28) 



G 



and when, a > 

~ z K(dz 



P 8 + \9 2 a 2 ' 



hence K{dz) = §§ e -( 2 P*)/ CT2 dz. 

The case of a = 0. If there is no Gaussian component, that is X t is 
drift process minus subordinator of bounded variation, then K(dz) = So(dz). 
The equation ^ allows also to identify e~ 9s ds J °° P(r+ > s)P x (r " < 
oo, -X T - G dz) using ([27]). (JTTJ) and (HU). Moreover, 

EX(1) Z.S'i , . 

p = P r - < oo = 1 = — = p . 29 

To identify ruin probability P x (r^ < oo) we need to find constant 
P(t^ < oo), which is given in the next theorem. 

Denote by p + (s) — P c (t _ < s) probability that the excursion above is 
shorter than s and let 

I* oo 

p(s, t) = / P(r+ +e < t) P € (t < s, -X T - G dz) 
Jo 

+P(T+<t)P e (T -<«,A T -=0) 

be probability that upper excursion above is shorter than s and the first 
consecutive excursion below 0, which is shifted downward by — e, is shorter then 
t. Note that p + (s) — p(s,oo). 

Theorem 2 (i) If X is a process of bounded variation, then 

f °° P(r+ > C)P(r - < oo, -X G dz) 



P(t c < oo) 
where 



1 - p + / °° P(r+ > C)P(r - < oo, -X T - G dz) ' 

e- 6s ds I P(r+ > s)P(t " < oo, -X.- G dz) 
Jo " 

= i/ OC ( 1_e "* ( ^) !T ^ (z)d " (30) 
(ii) If X is a process of unbounded variation, then 

P(t*<oo)= lim hm P+(&) - p(& ; C) . (31) 
V ^ fc^oo ,40 1 - p(6, C) 

Proof 'If X is a process of bounded variation then by [H Cor. VII. 5] is irregular 
for (—00,0). Since by drift assumption we excluded the case of compound 
Poisson process X, it means that is regular for (0, oo). Now point (i) follows 
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from Theorem Q] by taking x — and by using identities (j2U)l and (j3"5j) . In the 
proof of (ii) we adapt the idea from the proofs of main results of Dassios and 
Wu [7] and Dassios and Wu [5] . 

Let = 0. First we will define the sequence of stopping times: 

<7+ = inf{< > 5+ : X t < -e}, a~ = inf{t > S~ : X t = e}, 

S+ +1 = in£{t > a+ : X t = 0}, S' +1 = >a~:X t <0} 

and the processes: 

±= ( X t ±e if 6± < t < at 

Moreover, let 6 > and 

t 6 = inf{t > : t - sup{s < t : l ( x s <o)l(x t >o) = 1} > b}. 

Similarly, we define stopping times and t± for X^. Observe that: 

P(t£ < r+) < P(r C < n) < P(ri < t-). (32) 

We will find P(r^ < r^~) first by decomposing the path of the process X + into 
excursions above and below 0. Formally, probability of event Aj, that the first 
excursion below of length greater than £ is jth excursion and it happens before 
the first excursion above longer than b, equals: 

¥(A j )=p(bX) j - 1 (p + (b)-p(bX))- 
Summing over j = 1, 2, . . . gives: 

Similarly, 

P(r£ < T -) = P(r + > C) + P(T + < (34) 

Recall that process of bounded variation X was excluded from our considera- 
tions. Thus by regularity of for (-co, 0) it follows that is regular for (0, oo) 
(see [T71 Theorem 6.5, p. 142]). Straightforward consequence of this fact is that 
lirn e4 . P(T+ < C) = 1 an d that P(ri < r fc ~) and P(r^ < r fc +) have the same 
limits as e I 0. From (f3"2")l we then derive the assertion of the theorem. □ 

Remark 3 To find P(t c < oo) for the process (|26[) with a = one can use the 
following identity for z > (see [13l Corollary 4.5] and [TTJ p. 105 and Cor. 
7.5]): 

P(-X T - € dz,T ~ < oo) = p n 5 (dz) = i%(z)dz (35) 
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and then 



rOO rOO 

/ e~ es dsl P(r+ > s)Px(t " < oo,-X T - G dz) 
Jo Jo " 

-sft 1 -'"*)^!"*- 



Remark 4 The probability 1 — p(s,t) could be identified using double Laplace 
transform. Indeed, by jTTJ Ex. 6.7 (i), p. 176] and (flT]) . 



/3wjf J (l-p(8 i t))e-f >t e- u 'dtda = 1 - e -*<« e E e {e~* WX * 

&(8)X -, 

e ^^-l(-x ew >e) 
Ee*^- ' 

where e w is independent of X exponential random variable with intensity u>. 
Moreover, from 116 we know that 



$ (w) (j3 - u) 



(36) 



and that 

P(-X e G dz) = -^r^ (w) (dz) -wW^O^dz. (37) 

e <" $(w) 

Hence choosing function n(e) such that limit 

m(td) = lim (38) 

<40 n(e) 



exists and is finite, we have: 



l_p(M) 



oo />oo 



where 



= fl ,, 2 /^ M L, .u - (40) 



Then also we have: 



P(r< < oo) = lim 9(b,0-f,oo) 
y ' h^oo q(b,C) 

4 Cramer's estimate 

In this section we derive exponential asymptotics of Parisian ruin probability. 
Assume that there exists 7 > satisfying £2(7) = ip(— 7) = and that (p(9) 
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is finite in the neighborhood of 7. Then Ee 7X1 < 00 and we can define new 
measure P" 7 via ©. Define U^dx) = G\ 0) (dx) := e' x U(dx) and 

p OO 

H= / xU^idx), 



where u!f\dx) = / °° e-^ + ^F(H t € dx)cft for q > 0. Note from [3 that C/ 7 
is a renewal function of the ladder height process calculated on P~ 7 . Moreover, 
from [3] we have: 

gg z2 (o 1 /j) l gg(M) , r , 9 ^ 
- 9/3 = -^s -1 ^- (42) 

Note also that drift of X on P r is positive, since by (jHJ its Laplace exponent 
equals then 

(p^(9) = <p{-~t - 9) 
and hence ^'_ 7 (0+) = -^'(-7-) > 0. Thus 

P~ 7 (r+ < 00) = 1. (43) 

Theorem 3 We assume Cramer conditions, that is, that there exists a 7 € 
(0,oo) such that ^(7) = and that <p(6) is finite in the neighborhood 0/7. We 
also assume that support of H is not lattice when Tl(R) < 00. We have: 

lim e 7 *P x (r c < 00) = P(r c < 00)^^ + (1 - P(r< < oo))/ (c) (C),(44) 

jcfoo 7/i 

w/iere 

e- 9s /^( S ) ds = ^^ 



lf i6 (7 + <&(0)) 2 M 

and P(r^ < 00) is given in Theorem^ If fi = oo ; then LHS of {4$ * s understood 
to be 0. 

Proof By ([22]) and the Key Renewal theorem saying that U^(dx) on (0, 00) 
converges weakly as a measure to [i~ 1 dx (see [TTJ p. 188] and [3]) we have: 

lim e 72: P(^ < 00) = «(0, 0) lim e^ x U{x, 00) = Js(0, 0) lim UJx, 00) = ^2l22. 

rrfoo x\oo xfoo 7 

Moreover, by (gU and (g^ , 

e" 9s ds / P(t+ > s)P x (r " < oo, -X T - G dz) 
Jo " 

P(f+ < 00) 1 



-E 



-*(S)(X_ + -a;) 

e < 00 
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and by (|43l) . the Optional Stopping Theorem and [T71 Problem 5.5] we have: 



lim e~< x l 

a; "foo 



e T i ; ,t+ < oo 



= lim E~ 7 



e T ^ , -nj < oo 



= lim E~ 7 

X^OO 



-(*(9)+7)(X_ + -a!) 



= £_ 7 (O,$(0)+ 7 ) lim / e-^W^^-^Kvidy) 



K_ 7 (0, *(*) + 7 ) lim f e-W«-Hr)» fj i{x + dy) = ^-gg^L- 



□ 



( 7 + *(0))V 
This completes the proof in view of Theorem [T] 

5 Convolution-equivalent case 

We start from introducing class of Levy processes, which we will deal with in 
this section. Process X considered here have jump measure IT? in a class S" of 
convolution equivalent distributions, which could be formally defined as follows: 



Definition 1 (Class C^) For a parameter a > we say that distribution 
function G on [0, oo) with tail G = 1 — G belongs to class if 



(i) G{x) > for each x > 0, 
(ii) lim„. 



G -^" p = e Qa; for each and G is nonlattice, 

G(«) 



(m) limn^oo G ^ ra ^ = e Q if G is lattice (then assumed of span 1). 

Definition 2 (Class S^) We say that G belongs to class if 
(i) G G £ (q) for some a > 0, 
(ii) for some Mq < oo, we have 

lim £!M = 2M , 

u-^-oo G(u) 

where G* 2 (u) = 1 — G* 2 {u) and * denote convolution. 



(45) 



For all a € K such that the following integral is hnite we define the moment 
generating function 8 such that 



S a (G) = / e au G(du), 
Jo 
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where G is a finite distribution function. 

Recall that X = —X is a spectrally positive Levy process. Throughout this 
section we assume that for X we have: 



(0 



and 



(«) 



(Hi) 



% G 5 (q) for a > (46) 

f X Il^y) dy e (47) 
Jo 

£(a) < for a > 0; (48) 



e -9 5 a (fl")<l, where g = lira ^ ^ and S a (H) = £ (£Ti). (49) 



The first condition gives: 

% G 5 (q) for a > 0; (50) 

Latter condition (iii) has a force when a > 0; by drift assumption for a = this 
condition is automatically satisfied. 

We define relation: f(x) ~ g(x) iff lim^oo f(x)/g(x) = 1. 

Theorem 4 Under assumptions asymptotic Parisian ruin probability 

equals: 

P x (^ < oo) ~ EXx (nr C < oo) + (1 - P(r< < oo))/ (e) (C)) %(y) d», 

s f {e) (s)ds=^ (l - e-* (e)z ) B(z)dz (51) 
and 

S (*) = |^ + a jT e Q ^% (y) d„) . 

Above for a = i/ie quantity —(p(a)/a is understood in limiting sense and equals 
-(p'(0+)=EX 1 . 



where 

/•OO 



e 
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Proof From [TJJ Theorem 6.2] we have that: 



(r;- mi = /,-(().()i/-(,-. x! ^EX X ( 75 £L) / C,2) 



and that: 



lim P(X ? + — x G eb|r^ < oo) = B{z)dz, 



where 

The equations (j2"3")l and (|23|). dominated convergence theorem and Theorem [T] 
complete the proof. 

□ 



6 Examples 

6.1 General classic risk process (Op) 

For the process © we have %(dz) = XF(dz) and ip(9) = p9-\+\ e _ea F(dz 
By ill), 

P(r " < oo,-X T - G dz) = 

and hence 

/•OO 

e" es ds / P(r+ > s)P(r " < oo, -X r - G (fe) 

n JO 



(53) 



Probability P(r^ < oo) could be found using Theorem [2] with ([53]) and 

™, _ \ At/ 
p = P(r <oo = — , 
P 

where ^ = J^°yF(dy). Theorem [T] gives then Parisian ruin probability 
F x (t< < oo) for all x > 0. 

To find Cramer asymptotics, note that $(#) and 7 solve the equations: 

/•oo 

e-*( d >F(dz) = (A-p$(0)+0)/A, / e~< z F{dz) = (X+m)/X. 



>0 JO 

Moreover, from ([5]) we have: 

£(0,0) = v?'(0+) = p- Ai/. (54) 
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By (fl2|) using fact that ^(—7) = we have: 

(j, = A / ye 7V U^(y, oo) dy = X / ye iy F(y) dy 
Jo Jo 

and hence 

lim e^FJr^ < oo) = P(t c < oo) ( — [ ye lv F(y)dy\ 

xfoo \p - \v J J 

+(1-P(r< <oo))/< c )(C), 

where 

If F e for a > then P(r c < oo) is given in Theorem H 

6.2 Classic risk process ([1]) with exponential jumps 

Corollary 1 Assume that X t is Cramer- Lundberg risk process (QJ) with expo- 
nential claims F{dz) = £e~^ z dz, where £ = 1/v. Then 

to / C / ^ A -C 1 ^)* Z' 

whe 



^ < °°> - U- A(l-D) ' (55) 



D = 1 ~ f \ ^-e' {X+pi)t t' 1 h(2t^fpX)dt 
Jo V A 

cmd 7i (x) is modified Bessel junction of the first kind. 
Proof From (JTTJ) and (|TU)) (see also [TJ p. 63]) we have: 

k(0,0)U(x,oo) = Ae - ^)*. (56) 

Moreover, by lack-of-memory property of exponential distribution, distribution 
of undershoot of is also exponential with intensity £. Hence: 

[e-°° d s[nrt > s)F x (r - < oo,-X T - € *) = ^"(^-^^^ (57) 
and by |53|) 



/o 



e- e °dsj Q P(r+ > B )P(r - < oo, -X T - £ dz) — — > + ?) . (58) 
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Inverting the Laplace transforms (|57|) and f[55|) with respect to 9 (see [3]) gives: 

/•OO \ 

/ P(r+ > C)P«(7b" < ~ X r- e = —e- (E ^ )x D (59) 

Jo Pt, 

and 

j P(t+>C)P(t -<oo,-X t - edz) = -D. (60) 
Further, from Theorem [5] we have that 



A 



D 



P (- C < ») = fcf.w ( 61 ) 
+ p? 

Representation (fT6|) given in Theorem [1] and identities ([56l) . ([59]) and ([6lj) com- 
plete the proof of (f5"5j) . □ 
Parisian probability in this case was already identified in [TJ. 



6.3 Brownian motion with drift 
Corollary 2 Assume that 

X t = aB t +pt, 

where a, p > and i? t is a standard Brownian motion. Then 



,(t< <co)=e-< 2 ^' 2 ^ (62) 

$ ( £. / C | i p. Ol 

^lcrV 2 / cr V 2 

w/iere 

and M (.) is a cumulative distribution function for the standard normal distri- 
bution. 

Proof Note that: 

(r+ > 0P*(7-(T < - X r- G = °- 



Moreover, using 



VW^PP+^TT- (63) 



and (ITU1) we have 

W(x) = - f 1 - e^"^ 
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and hence from (|T7|) . 

P x (r - < oo) = e -( 2pCT " 2 ^. (64) 
The first passage probability T e + has inverse Gaussian distribution and 

Hrt < - V (2 ^ - ^) + e*"" V (-£ V? - ^ ) • (6« 
Similarly, 

P e (r - < 6) = e- 2p ™~ V f £ Vb - + Af(-?-Vb--^). (66) 

From (|65p ~ (|66p using continuity of Brownian paths note that we have: 
lim P + (b)- P (bX) = Um . P £ (r - < b)F(r+ > Q 



eio 1 - p(b, C) 1 - P e (r - < 6)P(r+ < C) 



^(fv / §)-fv / ^ 



(67) 



(68) 



(f>/§)+V?» (fx/?) 

where we use fact that W(— x) = ^(x). Note also that: 
hm^f^/fl = 4(2^Vl^(^)-^Vl +e_(f)2i 



6— >oo V 6 \ (T V 2 



Putting equations (|64|) and (|3Tj) into (|16|) of Theorem [T] completes the proof of 
(EH). □ 
Probability (f62j) was also given in [7]. 

One of the main goals of this paper was identifying Parisian ruin probability 
for more complex spectrally negative Levy processes. Such typical examples we 
analyze below. 

6.4 General classic risk process perturbed by independent 
Brownian motion 

We assume that: 

N t 



X t = x+pt-J2Ui + aB t 



i=i 
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where p, a > and B t is independent of Poisson process N t with intensity A and 
sequence of arriving claims {{7i} i=l! 2,.... Then ip{0) =p9-\+\ e - 9z F(dz) + 
^-6 2 and $(0) solves equation 



J e -*V>F(dz) = (\-p${0) - y$ 2 (6») + 9 



/A. 



Parisian ruin probability P x (t^ < oo) is given in Theorem[T]and all appearing 
there terms except P(r^ < oo) are given in Remark [2j 

Using Theorem [UJii) and Remark [4] we will find P(r^ < oo). Since X has 
a Gaussian component, then W^\0+) — and from Tauberian theorem (see 
e.g. [5J Theorem 1.7.1']) it follows that 

,. W^(e) 2 
hm = —z . 

eia e cH 

Thus n(e) = e and 

m(uj) 



ct 2 $(w)' 

Parisian ruin probability ¥(t^ < oo) is given in (|41|) for 



oo poo 



^""- ^"-'tM) - (70> 

Cramer asymptotics P^t^ < oo) is given in Theorem [3] where 7 solves 
equation: 

J\i'F(dz)=(\+py-^p) /A 
and /?(0, 0) = <^'(0+) = p — \u. Furthermore, 

H = A / ye™F(y) + a. 



Assume now that claim size has exponential distribution F(dz) = £e ^ z dz. 
Then by QJ] (see also [10]): 



2 



» x (t - < 00) = «(0, 0)tf(a;, 00) = 1 + *<T eiX , (71) 



i=l 
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where 



(i^ 2 +p) - v /(i^+p)2_2a2(^-A) 



(h& 2 +P) + v /(^a2 + p)2_2a2(^-A) 

-a 2 gf + 2pg 1 
<? 2 {qI - qI) ~ 2p(g 2 - Qi)' 



We have, 



/ e- 9s ds / P(r+ > s)P a (r - < oo, -X T - € cte) 
Jo Jo 



X. ) I 1 " 



(72) 

where $(0) solves equation 

2 

r 

2 " vw V" 2~ 



-* a (^) - ( £r +P ) ® 2 ( ) + ( A + 9 - PtMO) +£0 = (73) 



and by inverting this Laplace transform one can identify: J °° P(r+ > C)Pr ( T o < 
oo, —X- € rfz). 

We will find now 

P(r< < oo) = lim lim (l-p(h,C))-P ; (r ->6) 

for the exponential claim size with intensity £. Note that by (|T2"|) Laplace 

transform ^ ^1 — g^q^ J is then the only single Laplace transform that must 

be inverted at point s — £ to derive Parisian ruin probability P 2 (r^ < oo) 
(compare with ([70j) where double Laplace transform must be inverted). 
Note that 

P {b, C) = (1 - P e (r - < 6, X = 0))£ / P(r+ +e < C)^* <te 

u Jo 

+P(r+ <C)Pe(r - <6,X T - =0) 

and 

P £ (r - < 6,X r - - 0) - P e (r - < oo,X r - = 0) 

- / P z (r - < oo,X T - = 0)P £ (r- > 6,X h e <fe). (75) 
Jo 
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Thus from ([75]) and ([T7J Therorem 5.9, p. 122] we derive: 

l-P e (T ~ <b,X T - =0) /■« 
lim ^ = dw'(0+) - d / u(z) gi {b,dz), (76) 



where 



e (r > G dz) 



g x (6,dg) = lim y u ~ ' — -. (77) 

e|0 e 



If 

hm — — = o 2 , hm — — = 03(0) 

then from (ITHl 



to/ C ^ \ .92 +ff 4 (oo)- 93(00) 

P(T^ < OO) = ; r , 

92 + 34(00) 

where 

/•OO />00 

54 (6) = d(u'(0+) - / u(«)fli(6, dz))£ / P(r+ < C)e-« z dz 



and <7i(oo) = lim;,_j. 00 <^(&) (i = 1,3,4). The probability P(t+ < £) appearing 
in 54 could be found using Kendall's identity saying that if spectrally negative 
X t has has a density m(t, x) at x, then also has a density at s and 

P(r+ G ds) x , 

i = -m(«,x), (78) 

as s 

where in our case 



We will find functions <^ for i = 1,2, 3. 
Function gi(b, dz). Note that: 



\(r _ > 6, X b G dy) = P(rf- > 6, oB h - pb G dy) e - Ab 
+A / e~ At / P(rf- > t,aB t -pt G dz) 

/ P- z - ft +e(r - > 6 - t, X b _ 4 G dy) F(d/i) dt, 
Jo 



(79) 



where 

rf ~ = inf {t > : aB t - pt > e} 

and probability P 2 (t " > s,X s G dj/) appearing in (1791) could be derived using 
so-called ballot theorem: 

P z (t > s,X s e dy)/dy = m(s,y-z)-y -m(s - t,y)m(t,-z) dt; 

Jo s — t 

(80) 
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see [5]. Moreover, by the strong Markov property, 

z — e — p(t — s) 



P(rf - > t,aB t -pt G dz) = dz \ - J q p ( r f ~ € ds ) < /' 

where by Kendall's identity ([75| . 



y/t — sa 



Hence 

Hm P( T f->t > ^-p t 6d,) =B 
ej,o edz 



where 



Summarizing, 



^ / x 1 -P 2 * + 2zp /" 1 , , 

- (z,t) = —e ^ / —j=e 2^(ts) ds. 
27T in 



. gi (6,dy) = e - Ab S-(2;,6)dy 

+A / e - xt / H - (js,&) dz 



I P-,-j,(rJ- > 6 - i,X 6 _ t G dy) F(dfc) dt. 
Jo 

Function g3(b). Similarly, 
93(b) = e- Xb f E-(y,b)dy 

J — OC 

+A / e" At / E~(z,b)dz P- z - y (T Q > b-t)F(dy)dt. 



Quantity g2- Finally, 

P(r+ > C) = P(rf + > C)e~ AC 

+A / er xt / P(rf + > t, oB t + pt G dz) / P,_»(t+ > C - *) F(dy) di, 
J Jo 

(83) 

where 

rf + = inf {< > : aB t + pt > e] 
has the inverse Gaussian distribution: 

„, n, , . eds ( (sp—e) 2 

P(rf + G ds) = — exp > v y 
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lim *(r*+>t,*B t+P tedz) = 
e+o edz 



1 - p i t -2z P ft 1 



3 + (z,i) = — e I —=e *°Ht-s) ds 

27T In V,s 3 



Thus 



where 



and 

92 = e- x < [ E + (y,C)dy 

J — OO 

ft fO f oo 

+A / e~ xt 3+(z,C) / F z _ y (T+>(-t)F(dy)dt. (85) 

Jo J-oo JO 

The probability P z (T f > s) = P(t^ 2 > s) for z < appearing above could 
be obtained using Kendall's identity (175b . 

6.5 General classic risk process perturbed by a-stable mo- 
tion 

We will assume now that: 

Nt 

X t = x+pt-J2 u i + z t, (86) 
i=i 

where p, a > and Z t is a spectrally negative a-stable motion with a G (1, 2) 
independent of Poisson process Nt with intensity A and the sequence of arriving 
claims {f/i}i=i,2,...- Note that for a € (1,2) process Z is of unbounded varia- 
tion. Moreover, process X has no gaussian component and it does not creeps 
downward. Then ip(6) = pO — A + A f Q e~ 8z F(dz) + c9 a for some c > and 
$(0) solves equation 



e -*(9)« = (A - p$(0) - c$ a (#) + 0) /A. 

Then from |11) we have: 

oo 

P x (r - < oo) =k(0,Q)U(x,oo) = l-(l-po)£pJ (^(«+D* * M"*) (x), (87) 

n=0 

where Po = ^, M(dx) = ±i%) dy and K(dx) = £~ rd'/i)^ '"'' 1 " 

Finally, P x (t^ < oo) could be found using main representation given The- 
orem [TJ the equation (fT9|) and identity (|87|). Using Theorem [2jii) and Remark 
@]one can also identify P(r^ < oo) where similarly like in case of Brownian 
perturbation we have: (0+) = and 

(7W( f ) l 
hm — 
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Thus n(e) = e"" 1 and 

m(U!) = ; — ; ; -. 

v ; cT(a)$(w) 
Note also that for process (j8"6")l we have 

= + ( 88 ) 

Assuming that F(x) = o(x~ a ) for large x we derive that 11^ £ S^°\ From 
Theorem U we have 



t x (t c < oo) - (P(r c < oo)/EXi +P(r c = oo)/ (e) (C)) 



a(a — 1) 



-o+l 



as x — > oo, where the Laplace transform of f^(-) is given in (foTj) for defined 
in (1881). 
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